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Abstract 

We reconsider the problem of the stability of the thermohaline 
circulation as described by a two-dimensional Boussinesq model 
with mixed boundary conditions. We determine how the stability 
properties of the system depend on the intensity of the hydrologi- 
cal cycle. We define a two-dimensional parameters' space descrip- 
tive of the hydrology of the system and determine, by considering 
suitable quasi-static perturbations, a bounded region where mul- 
tiple equilibria of the system are realized. We then focus on how 
the response of the system to finite-amplitude surface freshwa- 
ter forcings depends on their rate of increase. We show that it 
is possible to define a robust separation between slow and fast 
regimes of forcing. Such separation is obtained by singling out 
an estimate of the critical growth rate for the anomalous forcing, 
which can be related to the characteristic advective time scale of 
the system. 
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1. Introduction 



The thermohaline circulation (THC) plays a major role i n the global circula- 



tion o 



1992; 



the oceans as pictured by the conveyor belt scheme ((Weaver and Hughe 



Stocker 



20011 ). The currently accepted picture is that the meridional 
overturning and the associated heat and freshwater transports are energeti- 
cally sustained by the action of winds and tides, controlling turbulent mix 



ing in the interior q 



Wunsch and Ferrari 



the ocean (Munk and Wunsch 



1998 : 



Rahmstori 



2003; 



2004). However, for climatic pourposes, many authors 



have succesfully assumed a dependence of the strength of th e THC from the 



Weaver and Hughes 


1992 


TziDerman et al. 


1991 


Marotzke 


1996 


Rahmstori 


1996; 


Gnanadesikan 


1999 


Stocker et al. 


2001 


)• 



The present day THC of the Atlantic Ocean is characterized by a strongly 
asymmetric structure. Deep convection is observed at high latitudes in the 
northern hemisphere. The water masses formed in the northern regions can 
be followed as they cross the equator a nd observed as they connect with the 



other major basins of the world ocean ((Weaver and Hughe 



2000, 


2002; 


Stocker et al. 


2001) 



1992; 



Rahmstori 



20011 1 . Idealized and realistic coupled GCM exper- 



iments have shown that such equatorial asymmetry may be a consequence 



equilibria ( 


Brvan 


198fi 


; Manabe and Stouffcr 


1988 


Stocker and Wright 


1991 




Manabe and Stouffei 


1999a 


Marotzke and Willebranc 


1991 


Hue-hes and Weaver 
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19941 ) . The equatorial asymmetry is also responsib 



the g l obal poleward 



1999; 



Stocker 



2nnnl 



leat transport 



Stocker et al.1 



(IBroecker 



1994 



e for a large portion of the 



Rahmstorf and Ganopolski 



20011 ). Consequently, large climatic shifts 



are often associated with important changes in the large scale oceanic cir- 
culation. On a paleoclimatic perspective, major cli matic shifts may be 



Bovle and Keigwin 


1987; 


Keigwin et al. 


19994 



Rahmstorl 



199.4 



1985 



2002). 



In fact, the THC is sensitive to changes in the climate since the North At 
lantic Deep Water (NADW) formation is affected by variations in air temper- 



ature and in precipitat ion in the Atlantic basin (jRahmstorf and Willebranc 



ma 



Rahmstorl 



1996). With respect to the present climate change, most 



GCMs have shown that the changes in radiative forcing caused by the on- 
going modification of the greenhouse gases in the atmosphere could imply 
a weakening of the THC. Large increases of the moisture flux and/or of 



1992: 



1999a 



els 



Vlanabe and Stouffer 



1993 



the northern Atlantic 


Weaver 


and Hughes 


Rahmstorl 


1997, 


1999b 




2000 


Wang et al. 



bJ). Moreover, models of di f ferent level of comp l exity, 



Tzipcrman and Gildor 



(Stocker and Schmittner 



1999; 



2002; 



1997 



Manabe and Stouffei 



Lucarini and Stone 



Schmittner and Stocker 



1999a 



b. 



rom box mod- 



2003a.jj), to E MICs 



1999|) to GCMs dStouffer and Manabe 



20001 ) have shown how the rate of in- 



crease of forci ng may be relevant for d e termin ing the response of the system. 



In particular 



Stocker and Schmittnerl ([19971 ) have performed a systematic 
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analysis of the stability of the meridional overturning circulation as a func- 
tion of the climate sensitivity and of the rate of CO2 increase. However they 
made no explicit reference as to the mechanism driving the response of their 
coupled model. 

In this work we study the THC stability usi ng a simplifi e d 2D Boussi- 
nesq ocean model, which has been presented in 



dime nsional models have been widely adopted (|Cessi and Youn g 1992; 



Art ale et al 



(2002) . Two 



Vellinga 



19961 ) and have 



proved their ability in de scribing the most relevant feedbacks 



of the system (jDiikstra and Neelin 



1999). Moreover, the low computational 



cost of such models permits extensive parametric studies. Our work wishes 
to bridge both in terms of methodology and results the studies performed 
with simplified models with the more physically sensible analyses performed 
with EMICs and GCMs. We explicitly analyze what is the role of the rates of 
changes of the hydrological forcing in determining the response of the system. 
In particular we determine, for a given initial state and a given rate of increase 
of the forcing, which are the thresholds in total change of the forcing beyond 
which destabilization of the THC occurs. The treatment of a wide range of 



temporal scales fo r the increase of the 



and continuously (jLucarini and Stone 



orcing a llows us to join on naturally 



2003al )bl) the analysis of quasi-static 



perturbations, which ha ye been usually addressed with the study of the bi- 



furcations of the system (IRa hmstorf 



Scott et al 



1999; 



Wang et al 



1999a; 



995 



1996 



Tit 



7- et a1 



the effects of very rapid perturbations 



Stone and Krasovskiv 



Rahmstorl 



1999; 



2002jjajh with the study of 



1996 



Scott et al. 



1999; 
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Wiebe and Weaver 



1999), which are usually differently framed. 
Our paper is organized as follows. In sectional we provide a description 
of the model we adopt in this study. In section |3] we explore the parameters 
space descriptive of the hydrology of the system by considering quasi-static 
perturbations. We determine which hydrological patterns are compatible 
with multiple equilibria and which hydrological patterns define a unique sta- 
tionary state. In section 0] we extend the analysis to time- dependent per- 
turbations. We analyze the temporal evolution of finite amplitude modifica- 
tions of the hydrological cycle that are able to destabilize the equilibria of 
this advective system. We propose a simple relation between an estimate of 
the critical rate of increase of forcing, which divides robustly slow from fast 
regimes, with an estimate of the characteristic advective time scale of the 
system. In section |S] we present our conclusions. In appendix A we present 
the dependence of the THC strength on the value of the vertical mixing 
coefficient. 



2. The model 

We consider the two-dimensional convection equations in the Boussinesq ap- 
proximation. The motion is forced by buoyancy gradients only: gravity is the 
only external force, while Earth rotation is not explicitly considered. Buoy- 
ancy gradients are generated in the model by imposing heat and freshwater 
fluxes at the top boundary. Such fluxes are assumed to be representative of 
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the interactions with the overlying atmosphere. 

We adopt a linearized equation of state for the sea water: 



p(T,S) = Po (l-aT + pS) (1) 

where the values of the coefficients of thermal expansion a = 8-10~ 4 °C~ l and 
haline contraction (3 = 1.5 ■ 10~ 3 psu~ 1 have been chosen in order to provide a 
good approximation over a quite large range of salinity and temperature. The 
linear approximation is commonly adopted in conceptual models. However, 
we note that analyses performed on simple models sho w that small nonlineari 



ties m ay induce self-sustained oscillations for the THC ((Rivin and Tziperman 



19971 ). Such nonlinearities in the equation of state might be especially rele- 
vant for the high latitude areas, since the well-known cabbeling effect occurs 
at low temperatures. 

The geometry of the model is descriptive of the Atlantic ocean, where 
we assume a depth of 5000m, an effective east-west extension of QOOOKm, 
and a north-south extension of 13600-K'm (120°) and a total volume V = 
4.08 x 10 17 m 3 . The only active boundary of the model is the air-sea interface. 

We select a relatively coarse uniform resolution with Nh x Ny = 64 x 
16 grid points, where Nh and Ny refer to the number of the horizontal 
and vertical grid points, respectively, in order to meet the computational 
requirements needed for performing a parametric study. 

A well known property of the THC system is the existence of regimes of 
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multiple equilibria. When equatorially symmetric surface forcing is applied at 
the surface, the equilibria of the system fall into three well-distinct classes. 
One class is characterized by the presence of two equatorially symmetric 
thermally direct cells, where the deep water is formed at high latitudes. 
Another class is characterized by equatorially symmetric salinity driven cells, 
where deep water is formed in the equatorial region. Equilibria belonging to 
these two classes are observed when either surface thermal or haline buoyancy 
forcings are largely dominant, respectively. When the two forcings have 
comparable intensity, multiple equilibria regimes - which constitute a third 
class - appear. In this case, the equilibria are characterized by the dominance 
of one overturning cell. If also geometry is symmetric with respect to the 
equator, the two equilibria have odd parity and map into each other by 
exchanging the sign of the latitude. 



a. Boundary conditions 

The boundary condition for the sea surface salinity is defined in terms of 
the imposed atmospheric freshwater flux F affecting the surface grid box of 
volume v — V/ (N H x N v ): 

d t S id=1 = — -Fj— , (2) 
v 

where we indicate the value of the bulk variable S at the grid point with 
Sij and the value of the interface variable F at the grid point (z, 1) with F{. 
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We underline that, in order to simplify the expressions, in equation J2J (and 
in the following ones) we have not explicitly adopted a discrete notation for 
the time variable. 

ected the contribution 



19961 ) 



We emphasize that in expression (J2J) we have neg i 
in terms of mass of the freshwater flux to the ocean (|Marotzkd 

We divide the water basin into three distinct regions by using a suitable 
analytical expression of the freshwater flux. The equatorial region (region 
E) is characterized by a net atmospheric export of freshwater, while the 
northern and southern high latitude regions (regions N and S, respectively) 
are characterized by a positive atmospheric freshwater budget. We then 
consider the following functional form for the surface freshwater flux: 



Ft = -$i cos (2m/N H ) 

7T 



(3) 



where: 



$ 5 , i < 1/4 N H 
$tv, *>3/4N H 

§ E = -l/2($jv + $ 5 ), l/4JVff<i<3/4JV fl 



(4) 



The definition of Fi is such that $jv arid $5 are respectively the value of 
the total net atmospheric freshwater fluxes into regions N and S, while 
is constrained in order to have conservation of the salinity of the ocean. 
This latter condition is needed to allow the system to reach equilibrium 
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states. Therefore, the ocean conserves its average salinity So = 35 psu. Since 
the atmospheric freshwater budgets fo r the Atlantic ocean are though t to 



19751), we 



be positive for the regions N and S (jBaumeartner and Reichell 
consider the case $at, $5 > 0. Three relevant examples of surface freshwater 
flux are depicted in figure [T] 

The sea surface temperature is restored to a time-independent climato- 
logical temperature field Tj with a newtonian relaxation law: 



dtTij = i — X \Ti — Tj J= i] 



(5) 



where the constant A describes th e efficiency of the process. Such very simpli- 



fied ocean-atmosphere coupling ([Marotzke and Stone 



1995 



Marotzke 



19961 ) 



synthetically describes the combined effects of radiative heating-cooling and 
of the atmospheric latent and sensible heat meridional transport. The clima- 
tological temperature profile T; profile is set with the following equatorially 
symmetric analytical form: 



Ti = T - AT cos (2%i/N H ) 



(6) 



where 2AT is the imposed equator-to-pole temperature gradient and To is 
the average value of Tj. In accordance to definition of the E, N, and S 
regions with respect to their atmospheric freshwater budget, we observe that 
regions N and S are cold at surface, i.e. Ti is lower than T , while the region 
E is warm at surface, i.e. Ti is larger than T . 
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Since in our work we explore the stability of the THC with respect to 
the hydrology of the system, we keep fixed the parameters determining the 
restoring temperature profile Tj. We set T = 15°C, since it represents a 
reasonable average surface climatological temperature, and we choose AT ~ 
23.5°C, which corresponds to forcing the surface equator-to-pole temperature 
gradient to be ~ 47°C. Such a choice also implies that the average T is 
~ 0°C in the two high latitude regions N and S and is ~ 30°C in the 
eq uatorial re g ion E . Furthermore, following the parameterization proposed 



bv iMarotzkd ( 1996 ). we choose the restoring constant A ~ (ly)^ 1 , which is 
reasonable for the description of the temperature relaxation of the uppermost 
5000/16 m ~ 300 m of the ocean. 



3. Quasi-static hysteresis and multiple equi- 
libria 

a. Symmetric case: &s = $jv 

We start by considering the existence of a region of multiple equilibria in 
the space of parameters that defines the hydrology of the system, i. e. the 
($s,$at) plane. The surface thermal forcing is kept fixed. In the space of 
parameters that we are considering, symmetric forcings constitute the bisec- 
trix of the ($5,$jv) plane. The allowed circulation patterns change when 
adjusting the parameters along the bisectrix, as discussed in section |21 For 
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a weak hydrological cycle ($at = &s < ^m/)> we observe a stable symmetric 
circulation with downwelling at high latitudes. If the hydrological cycle is 
very strong ($^r = $5 > $ sup ), we obtain a stable symmetric circulation 
with downwelling of warm, saline water at the equator. In the intermediate 
regimes ($m/ < &n = &s < &su P ), the system has multiple equilibria. For 
later convenience, we define $ a „ = 0.5 + & SU p), which corresponds to an 
average hydrological cycle. The described equilibria are depicted in figures 
Et), 12b), Ek), arid Eh") respectively. The two points ($5 = $ in /,$Ar = $; n /) 
and ($5 = $ sup , $tv = are bifurcation p oints in the one-dimensional 



subspace = $5 (jDiikstra and Neelin 



1999T ) . We consider the northern 



sinking equilibrium state (note that northern and southern sinking patterns 
are equivalent in these terms) realized for $5 = $tv = &av a s the reference 
state of the system. 



b. General case: $s 7^ $jv 

The study of the multiple equilibria states can be extended to the case of 
non-symmetric forcings, i.e. $n. We wish to obtain an estimate of the 

shape of the domain T in the ($5, $at) plane where the system has multiple 
equilibria. 

Since, apart from the freshwater flux boundary conditions, the model is 
wholly symmetric with respect to the equator, we expect that any property 
of the system is invariant for exchange of $a? and $5, so that we have that 
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r is a priori symmetric with respect to the bisectrix. 

A fundamental property of the region T is that, if we start from a point 
belonging to T and change quasi-statically $5 and $tv along a closed path so 
that the point ($5, <3?tv) remains inside T, we get back to the initial equilib- 
rium state. Instead, if the closed path crosses the boundary of T, the initial 
equilibrium state may not be recovered at the end of the loop, since the final 
state depends on the path. 

As a starting position of the loop, we consider a northern sinking equi- 
librium corresponding to one of the stable states of the point &i n f < $s = 
&n < &sup on the bisectrix. By definition, such point belongs to T. How- 
ever, the initial position of the loop - if belonging to T - is not relevant in 
determining the shape of the region of multiple equilibria. 

We increase the value of $at at a given slow constant rate r s over a time 
to and then decrease it back to the initial value at the rate — r s . By slow 
rate we mean that r s <C <& a v/T] we select r s = $ a „/100r. If the initial state 
is not recovered at the end of the integration, we deduce that the path has 
crossed the boundary of T. By bisection, we can determine the critical value 
£q™*, which determines (^s, $tv + r s • = &n + A$^ J< ) as belonging to the 
boundary of T. A schematization of this procedure is depicted in figure 0^). 
By changing the initial point along the considered segment of the bisectrix, 
we are able to define the whole boundary of T above the bisectrix with a good 
degree of precision. Then, the symmetry properties of T allow us to easily 
deduce the portion of its boundary lying below the bisectrix. In figure El we 
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present the estimate for the the boundary of the bistable region T obtained 
following this strategy. 

We conclude that the boundary of T is constituted by the bifurcation 
points of the system. We note that in the more general case of the 2D 
$ N ) plane, ( $,<? = $ <WLf ,$y = $ in/ ) and ($ s = $ sup , § N = $ sup ) result 



to be cusp points ((Diikstra 



200l|). 



4. Effects of transient perturbations 

The analysis we have performed captures the equilibrium properties of the 
system, but is not sufficient to gain insight on the response of the system 
to transient changes in the forcings, which in general can range from instan- 
taneous to quasi-static perturbatio ns. Such sort of problems has been first 



investigated in the seminal paper by 



Stocker and Schmittnerl f)1997l ) with the 



purpose of determining how efficiently the negative feedbacks of the system 
can counteract external perturbations, depending on the temporal patters of 
the destabilizing forcings. 

In our case, it is reasonable to expect that, if we change <3>at at a finite 
rate, the system can destabilize before reaching the boundary of T. In fact, 
a fast perturbation can overcome the ability of the advective feedback to 
stabilize the system. We also expect that the faster the perturbation, the 
stronger such effect, i.e. the smaller the total perturbation required to obtain 
destabilization. 
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As in the previous case, the analysis starts by considering as initial equi- 
libria the northern sinking states under symmetric forcing. In this case, we 
increase the value of & t a constant rate 77 over a time to, we let the 
system adjust for a time lOr, so that transients can die out, and then de- 
crease $at back to the initial value at the slow rate — r s , corresponding to a 
quasi-static change. We schematically depict such strategy for three differ- 
ent values of 77 in figure |3Jd). If at the end of the process the initial state is 
not recovered and southern sinking state is instead realized, the system has 
made a transition to the other branch of the multiple equilibria. Depending 
on the choice of 77, we obtain different values for the critical perturbation 
causing such transition. Along the lines of the quasi-static analysis, varying 
the initial point, we can obtain for each value of 77 a curve which describes 
the critical perturbations. 

In figure El we report the curves obtained by selecting, from fast to slow, 
77 = 00, 77 = $ au /r, 77 = $ai;/3r, 77 = $ at) /10r, and 77 = r s = $ a „/100r. 

If we select 77 = r s , we obtain by definition the previously described 
upper branch of the boundary of T depicted in figure since the presence of 
the relaxing time lOr is not relevant. 

On the other extreme, if we apply instantaneous changes of $jv ( r / — 00), 
we obtain information on the minimum change in $jv that is needed to 
destabilize the system for any initial state having symmetric surface forcing. 
In fact, the corresponding curve is the closest to the bisectrix. 

Considering intermediate values of 77, we obtain consistently that the 
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curves of the critical perturbations lye within the two extremes obtained 
with rf = r s and 77 = oo. Moreover, we have that the curves are properly 
ordered with respect to the value of 77, i.e. the smaller 77, the closer the 
corresponding curve to the upper branch of the boundary of T. 

Previous studies, albeit performed with coupled EMICs, obtain a qualita- 
tively similar dependence of thresholds on the rate of increas e of the forcings 



(Stocker and Schmittner 



1997; 



Schmittner and Stocker 



1999), while in other 



studies (on GCMs) where the full collapse of the THC is not obtained, it is 
nevertheless observed that the highe r the rate of increase of the forcing, the 



larger the decrease of THC realized (jStouffer and Manabe 



1999). 



The most important result is that we can identify two separate regimes. 
If the surface forcing changes with a rate faster than $ au /r, the response of 
the system is virtually identical to the case of instantaneous changes. On 
the other side, if the rate of change is smaller that $ a „/10r, the response 
of the system is very close to the case of quasi-static perturbations. The 
curve corresponding to 77 = $ a ^/3r is geometrically about midway between 
these two regimes and is not patched to either. We have that the response of 
the system to varying external perturbations dramatically changes when the 
time scale of the variation of the external forcing changes by only one order 
of magnitude. Therefore, we can interpret r c = $ a „/3r as an estimate of the 
critical rate of change of the hydrology of the system. 

It follows that with r c we identify a relation between the class of changes 
in the external forcing that distinctively affect the stability of the system and 
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the internal time scale of the system. 

The previous results prove to be robust with respect to changes in the 
freshwater forcing of the initial states. As an example, in figure[7|we show the 
results of a similar analysis referring to initial states having non-symmetric 
surface freshwater forcing ($tv = 2$s). 

We explored the behavior of the system in the ($5, $at) plane considering 
only changes in $jy for computational convenience. Nevertheless, coherent 
results can be obtained by changing both parameters. In this case rj has to 
be interpreted as the sum of the absolute values of the rates of change of the 
two parameters. 

Changing the parameter AT implies a change in the values of <&i n f and 
$ snp , since the meridional gradient of the total buoyancy forcing is changed. 
However, the response of the THC to transient perturbations does not change 
qualitatively. Instead, figures EE must be rescaled linearly with the proper 
values of and $j n /- Such linear relation is a direct consequence of the 
use of a linearized equation of state for the sea water. 

5. Conclusions 

This work provides a complete analysis of the stability of the ocean system 
under examination with respect to perturbation to the hydrological cycle. 
We have provided a simple description of the profile of the net freshwater 
flux into the ocean which is fully specified when only two parameters, which 
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are related to the total freshwater budget of the two high-latitudes regions, 
are specified. 

We have first analyzed the bifurcations of the symmetric system, which 
might be taken as the prototype of a system that has equal probabilities of 
falling in two different equilibrium configurations. We have found that the 
system is characterized by two bifurcations, which delimitate a domain of 
multiple equilibria. 

We have then extended the study to the general case where asymmetries 
in the hydrology are considered. We have produced a two-dimensional sta- 
bility graph and have pointed out the presence of a region V where multiple 
equilibria are realized. Our results summarize the information that can be 
obtained with multiple hysteresis studies. 

In this study we emphasize that the rate at which a perturbation in 
the hydrological cycle is applied to a simple model of the THC may dra- 
matically affect its stability. When general time-depends perturbation to 
the hydrological cycle are considered, we obtain that the shorter the time 
scale of the forcing, the smaller the total perturbation required to disrupt 
the initial pattern of the circulation. The observed relevance of the tempo- 
ral scale of the forcing in determining the response of our system to per- 



turbations affecting t 


e st 


abi 


itv of the THC agrees wit! 


i the findings of 
) for box models, 


TziDerman and Gildor 


(20(1 


2). 


II 

Lucarini and Stone (2003allb 


of Stocker and Schmittner 


1997) and Schmittner and S 


tocker ( 199S ) in the 


context of EMICSs, and of ( 


Vlanabe and Stoufferil999a). 


Manabe and Stouffer 
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(1999b) 



Manabe and Stoufferl ([2000D . and 



Stouffer and Manabd f)1999f ) in 



the context of GCMs. Moreover, the saturation and patching effect observed 
for slowly and rapidly increasing perturbations, which allows the definition 



of slow and fast regimes, resp ectively, agrees qualitatively with the findings 

3) 



of 



Lucarini and Stone 



2003a b) for box models, and resembles some of the 



results - albei t obtained with a coupled mode l and c onsidering CO? increases 



prese nted in 



Stocker and Schmittnerl (|1997l ) and in 



Schmittner and Stocker 



(J1999J). 

The main conceptual improvement we propose in this work is the exis- 
tence of a relation between the critical rate of change of the forcing and the 
characteristic advective time scale of the system. 

We notice that the advective time scale results to depend on the inverse 
of the square root of Ky, as discussed in appendix A. Therefore, a very 
important consequence of this analysis is that the efficiency of the vertical 
mixing might be also one of the key factors determining the response of 
the THC system to transient changes in the surface forcings. Future work 
should specifically address the details of the functional dependence of the 
critical rate on the the vertical diffusivity. 

Other relevant improvements to the present study could be the adoption 
of a more complex ocean model, descriptive of other ocean basins, as well 
as the consideration of a simplified coupled atmosphere-ocean model, where 
the effects of the radiative forcing can be more properly represented. 
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Appendix A. Relevance of the vertical diffu- 



sivity K 



v 



The vertical diffusivity Ky, or, equivalently, the diapycnal diffusivity K D , is 
the criti cal parameter controlling the maximum THC strength ty ma x m ocean 



models ( Brvan 



1987 



Wright and Stocker 



1992). On the other h and, an esti 



mate of its value in the real ocean is a subject of current research (jGregg et al 



2003|). Scaling theories proposing a balance between vertical diffusion and ad- 



vection processes suggest, in the case of three-dimensional hemis pheric model 



of th e Atlantic ocean, a power law dependence \I/ r 



1999; 



Dalan et al. 



Ky 3 (|Zhang et al 



20041 ) . In the cas e of two-dimension al models, the ex 



pected dependence is \l/ r 



Ky 2 ( 



20001) . This relation is 



verified in our model in the range 0.6cm 2 s -1 < Ky < 4.0cm 2 s' 1 (figure IHJ). 
In this study, the value of Ky has been selected so that the corresponding 
northern sinking equilibrium state characterized by an hydrological cycle de- 
termined by $5 = $tv = &av has an overturning circulation ~ 30Sv . With 
this choice of Ky we can define a characteristic time scale r for the sys- 
tem as r = V/^ m ax ~ 350y. Similarly, if we change Ky over the range 
shown in figure |H1 the advective time-scale would range over the interval 
150y < t < AOOy. These considerations will be useful in the final discussion 
of our results. Given the parameters chosen for our simulations, our model 
integrations estimate ~ 0.04S'f and ~ 0.73Sv, so that <3> a „ results 
to be ~ 0.39Sv. 
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Figure 2: Circulation patterns obtained in the non-bistable region for sym- 
metric forcing. Transport is expressed in Sv. a): _s = _jv = 0.1_- ai) = 
0.5_. n/ . b) : _ s = _ jv = 2_ at) = 1.1 _ flU -. 
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Figure 3: Circulation patterns obtained in the bistable region for symmetric 
forcing with $5 = $tv = &av Transport is expressed in Sv. a) Northern 
sinking pattern, b) Southern sinking patterns. 
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(a) Schematization of the <£>at time-dependence for the 
three cases of loop within T, loop crossing the border of T, 
and loop reaching the border of T. 



Rapidly increasing forcing 

More rapidly increasing forcing 

— Instantaneous forcing 




Time 

(b) Schematization of the $jv time-dependence for three 
cases of transient forcings. 



Figure 4: Quasi-static and transient changes in the value of the parameter 
<5>n 
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Figure 5: Stability graph of the system in the space ($5, $at). The thick 
black line delimitates the bistability region T. Along the diagonal, solid lines 
represent the bistable states having antisymmetric circulation patters, while 
the dashed line represent represent the stable symmetric circulation patters. 
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Figure 6: Critical forcings for the collapse of the northern sinking pattern in 
the space ($5, $at). Various temporal patterns of the forcings are considered. 



35 



- 


1 


i 
























<^ 

r ¥ 
* > 












✓ 

✓ 

✓ 

✓ 






* X 


► X 








✓ 

✓ 

✓ 

✓ 








1 f 

I / 4 

1 *i 
■ ^ # 










✓ 

✓ 

✓ 










\ * J 

1 

IV 










r=r ~ 1 $ /100 x 

f s av 

- - ■ r ~ 1 <I> /1 x 

f av 

r -1$ /3 x 

f av 

- - ■ r ~ 1 <£ /1 x 

f av 

Inst, flux changes 

O * inf 
w sup 

<£>=<!> - mult. eq. 

N S ^ 

$ N =$ s - sym. circ. 

▲ Bifur. <I> =2* 

▼ N S 

In. states <I> =2* 

N S 








✓ 

✓ 

✓ 

✓ 








✓ 

✓ 

✓ 








4 

* 

* 

* 

* 








✓ 

✓ 

* 












i 











0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

<!>/<!> 

S av 

Figure 7: Critical forcings for the collapse of the northern sinking pattern in 
the space ($5, $at). Various temporal patterns of the forcings are considered. 
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Figure 8: Dependence of maximum value of the THC on the value of the 
vertical diffusivity for T = 15.0°C, AT = 23.5°C, and $jv = $5 = 
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